Glueball filled with the quark field as a model of nucleon 
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Abstract 

On the basis of the non-perturbative Heisenberg's quantization scheme 
and using some simplifications and assumptions the reduction from the 
gluon-quark Lagrangian to a scalar-fermion Lagrangian is made. The cor- 
responding field equations have a regular spherically symmetric solution 
which can be interpreted as a simplified model of nucleon. Some prop- 
erties of such a "nucleon" are discussed: the presence of a mass gap and 
spin ^ h/2. The mass gap describes the difference between the "nucleon" 
and the "nucleon" without quark field (a scalar glueball). The presence 
of spin ^ h,/2 can be connected with the spin problem of nucleon. 



1 Introduction 

One can assume that removing the quarks from a nucleon gives us a glueball. 
In this case one of the methods to prove a glueball model is by filling up it 
with fermions (quarks) and comparing it with nucleon properties. In Ref. PP a 
scalar model of glueball is presented. This model is based on a non-perturbative 
quantization technique which was applied by Heisenberg for a non-linear spinor 
field. In our approach we have applied this idea for quantum chromodynamics. 
The essence of this non-perturbative technique is that in the first approximation 
the field degrees of freedom {A^ gauge potential, B = 1,2, ...8 is the color 
indices, ji = 0,1,2,3 is the world indices) can be reduced to scalar fields 4)^. 
The components (j)°-{a = 1,2,3) and (jf^ (m = 4,5,6,7,8) of the scalar field 
have a different dynamical behavior. The numerical analysis shows that the 
corresponding field equations for 0"'™ have regular solutions only for some values 
of parameters of the model. The solution looks like a bag, into which we would 
like to put quarks (fermions). 



2 Short description of the model 

In this section we would like to derive the field equations describing in an approx- 
imate manner quantized fermion and SU(3) gauge fields. Strictly speaking the 
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Heisenberg's quantization procedure |2] is based on an infinite set of equations 
relating all Green's functions which are similar to Dayson-Schwinger equations 
but on the non-perturbative language. Of course such a set of equations is math- 
ematically extremely hard for solving and we need for some physical reasonings 
for a closure. In the approach presented here we consider only 2 and 4-point 
Green's functions. In this case we can average the quantum chromodynamics 
Lagrangian and derive equations describing 2 and 4-point Green's functions. 
For this purpose we assume that a 4-point Green's function can be presented 
as some sum of the products of 2-point Green's functions, and 2-point Green's 
function can be decomposed as the product of scalar fields. The Lagrangian of 
the SU(3) gauge field interacting with quarks is 



(L C^Q ^" (La 



4 z^'' 



(1) 



where Cg q are the gluon and quark Lagrangians correspondingly; F^^ — ly^j 
d^A^ + Gf^^^A^Aj^ is the field strength operator; C, £» = 1, . . . , 8 are 
the SU(3) color indices; G is the coupling constant; f^^^ are the structure 
constants for the SU{3) gauge group; A^ is the gauge potential operator; A^^ = 

L^A^; are the generators of the SU(3) gauge group; -0 is the quark field 
and M is the mass matrix of quarks. In order to derive equations describing 
the quantized field we average this Lagrangian over a quantum state | Q) 



Q 



Q) = (£) = (£ 



(2) 



where 



- ic/^^^ (.) - d.A^^ix)) A^^{x)A^'^{x)) 

(d^ + iGA^) $ 



(3) 
(4) 



At first we will calculate i^^g^ following Refs. One can see that schemat- 

ically we have the following 2, 3, and 4-point Green's functions: (^(dA)"^^, 

{{dA)A^) and ((^)*)- We suppose that the odd Green's functions can be 
written as the following product: 



A^{x)A'^{y)Ai^{z)^ w (^A^{x)A'^{y)j (^Ai^iz)j + (other permutations) = 

(5) 
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as we suppose that (Af (x)) = 0. Further we suppose that a 2-point Green's 
function can be presented in the so-called one-function approximation IB1 as 



(6) 



where (j)^{x) is the scalar field which describes the 2-point Green's function. 
These two assumptions are similar to the quantum harmonic oscillator where 
(a;) = but (x^) ^ 0. The 4-point Green's function can be written in a one- 
function approximation as in the symmetrized product of corresponding two 
2-point Green's functions 



(7) 



A^{x)A^{y)Ai;{z)Afiu)) « {A^{x)A'^ (y)) {A^{z)Af{u)^ 
A^{x)A^iz)) (A^,{y)Af{u)) + (A^ix)Af{u)) (l^(y)l^(z; 



Taking into account these expressions for the 2,3, and 4-point Green's functions 
we can derive the effective Lagrangian £0 
describes 2 and 4-point Green's functions 



Cg ) for the scalar field (j)^ which 
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where the indices a — 1,2,3 are SU(2) indices and m — 4,5,6,7,8 are the 
coset SU(3)/SU(2) indices; are some constants. We can add the term 
— A2(</!'™(/'™)^/4, which is inessential for the dynamics (but essential for the 
finiteness of the energy) , and the Lagrangian becomes 
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The next step is calculation of {Cq 



(10) 



We see that the fermion (quark) field ip does not have any strong self-interaction. 
Therefore we can suppose that in the first approximation the dynamical behavior 
of the quantized quark field ip is similar to the dynamics of classical field 



(11) 



Here, we have to note that according to our previous assumptions (^^) = and 
consequently the second term in Eq. Hll|l is zero but the interaction between 
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gluon and quark fields must exist. To introduce the interaction between (f>^ and 
tjj we have to do some physical assumptions about this interaction and insert 
this term into Lagrangian by hand. For example, one can assume that the quark 
field interacts with fluctuations of the gluon field, 



(12) 



where for the redefinition of a' we have used Eq. (jSJ . This assumption is equiv- 
alent to cutting off the infinite set of equations connecting all Green's functions 
for exactly quantized gluon and quark fields. Thus the quark Lagrangian is 



(13) 



Hereafter, we omit flavor indices to simplify notation. Finally we have the 
following effective Lagrangian describing the non-perturbative quantized qluon 
and quark fields interacting amongst themselves 
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Ai 



20i 



v^o 



(14) 



M - «( 



with the following field equations 



20 + Al (0 - 00 00 ) 



[i-f^df, - (M - a0^0^)] iIj = 0. 



(15) 

(16) 
(17) 



3 Initial equations 



In this section we would like to present the solution which later we will view as a 
simplified model of the nucleon. We search solution in the following spherically 
symmetric form; 



0-(r) 
ip{t,r,0,(p) 



0(r) 
fir) 



-iEt 



a = 1,2,3; 

m = 4,5,6,7,8; 



g{r) COS0 
\ g{r) smOe-^f J 



(18) 
(19) 

(20) 
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where r, 0, ip are spherical coordinates. Once again we would Hke to remind that 
in our approach we suppose that the components of scalar fields with different 
indices a — 1,2,3 and 0™, m = 4, 5, 6, 7 have different dynamical behavior. 
After s ubstituting Eqs. |(TSJl-l|2ni) into Eqs. ^^-^7^ gives 



+ 



/" + -/' 

r 

h' + -h 
r 

9' 



f 



{h'-9') 
ih'-9') 
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-h 



(21) 
(22) 

(23) 
(24) 



where 20g0o — m? and 2(/)™(/)™ — fi^; to, fi are some constants which will be cal- 
culated by solving Eqs. |[^ - (P^ : and constants Ai,2 are redefined Ai,2/2 Xi_2- 
We redefine 0(r)/0(O) ^ (j}{x), /(r)/0(O) ^ f (x), m/(j}{<S) fi/cjjiO) 
M, ft(r)/03/2(O) /i(x), g{r)/(b'/^0) -> g(x), a/(/.(0) ^ a, E/m ^ i?, 
M/(j){0) M and introduce the dimensionless radius x — r0(O). After this we 
have the following set of equations: 



0" + -0' 

X 

f" + -r 

X 

h' + -h 

X 



— TO — 



ih'-9') 



= f 



5p + M-|(/2 + 02) 
-h[E-M + ^{f + cj^^) 



(25) 
(26) 

(27) 
(28) 



where d{- ■ ■ )/dx = (•••)'• Evidently these equations are too complicated to be 
solved analytically. Preliminary numerical investigations show us that the set 
H25l) - (|28|l does not have regular solutions for arbitrary choice of the parameters 
TO, fi, and E. We will solve equations H25l) - H28l) as a nonlinear eigenvalue problem 
for eigenstates (/)(a;), /(x), h{x),g{x) and eigenvalues to, /j,, E, i.e. we have to find 
parameters m,fi,E to provide existence of regular functions 4>{r), f{r), h{r) 
and g{r). We will search for the regular solution under the following boundary 
conditions 



0(0) = 1, 0(oo) = to; 

/(O) = /o, /(oo) = 0; 

h{0) = 0; 

5(0) = go- 



(29) 
(30) 
(31) 
(32) 



The boundary conditions for (t>{x) and f{x) allow us to say that these functions 
resemble interacting kink and soliton. 
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For the regular solution which will be presented below we assign the following 
values: 0(0) = 1.0; /(O) = VoI;g{0) = 0.1; G = 2.0; M = 0.2; Ai = 1.0; As = 
0.1; a = 1.0. 



4 Numerical solution 

In accordance Ref. [I], we use the following numerical method for solving Eqs. (^5] 
l|28(l : we take a zero approximation for the functions f{x), h{x), g(x) (which are 
fo{x) — 0.6/ cosh^(x/4) and ha{x) = go{x) — 0) and solve equation in the 
following form: 

0'/ + -01 = 01 [/2 + A2 (0? ~ ml)] (33) 



where mi is the first approximation for the parameter m, the boundary con- 
ditions are H29|) and the function 0i (x) is the first approximation for the func- 
tion 0(a;); 0i(a;) and mi are the eigenfunction and the eigenvalue for Eq. H25|l 
correspondingly, with the boundary conditions (|29|) . Obtaining the regular 
solution 01 (x) we can substitute it into Eq. H26() with zero approximations 
ho{x) = goix) = and solve the equation 



fi'^-fi 



fi [0? + Ai (/? - fil)] 



(34) 



with the boundary conditions (|30() . fi{x) and /ii are correspondingly the eigen- 
function and the eigenvalue for Ea. (P()|l . with the boundary conditions 

Further having the first approximations 0i(x), fi{x) we can substitute them 
into Eqs. ^ and ^ 



X 



= .91 



i? + Af-^(/f + 0?) 



-hi 



E-M+^{f,+c^l) 



(35) 
(36) 



where hi{x), gi{x), Ei are eigenfunctions and eigenvalue for Dirac equations H35|l 
and H36I) with the potential ^(/i + 0i). Thus having the first approximations 
for (j)i{x), fi{x), hi{x), gi{x),mi, ^i, El one can substitute these quantities into 
Eq. (ESJ, 



- mi - 



{hi - gl) 



(37) 



This equation gives us the second approximation 02 (x) for the function (j){x), 
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and so on. Thus, on the r step we have 



2 

4>'l + -(t^'i 

X 

K + Ik 



9i 



fli + Ai 



(hli-gli) 



(38) 
(39) 

(40) 
(41) 



where the functions fi-i{x), hi^i{x), gi-i{x) are defined in the preceding step, 



// 2 



E^-l ~M + - {fl, + 0ti) 



Such procedure has to lead to precise eigenfunctions 

4>:{x)^4>*{x), f*{x)^rix), h*ix)^h*ix), g*{x) 
with precise eigenvalues 

m* m*, ^l* m*, E* E* 



(42) 

(43) 
(44) 



g*ix) (45) 



(46) 



where (/)* , f* , h* , g* are the regular solutions of Eqs. which can be 

regular only for the values of parameters m* , fi* , E* . For the other values of 
parameters nii ^ m* , ^ IJi-*,Ei ^ E* the functions 4>i ^ , fi ^ f* , hi ^ 
h*,gi g* are singular and imply infinite energy. 

4.1 More detailed description of the numerical calcula- 
tions for every step 

In this section we would like to describe more carefully the procedure of the 
numerical solution of each equation of the system H38() - (|44(l . 

At first we would like to describe more carefully the numerical solution of 
equation H33|) . For this purpose, we choose the zero approximation fo{x) as 
follows: 

foix) = (47) 
cosh I 

Typical solutions for the arbitrary value of nii are presented in Fig. ^ We see 
that when mi < m* (ml is an eigenvalue which gives us an eigenfunction (jiKx)) 
the solution (pi {x) is singular and closer to the singularity the equation has the 
form 
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10- 



8- 




Figure 1: The typical singular solutions (/)(a;) of Eq. The solution is pre- 

sented for i = 4 step. 



< « X2<I>1 (48) 

Consequently the solution is 

(49) 

V M xq- X 

where xq is some constant depending on mi. On the other hand, for mi > ml 
the solution is also presented in Fig. ^ and the corresponding asymptotical 
equation is 

(f>'l{x) + ^cj)[ « - (X^m^) (50) 
which has the following solution: 

sin ix\/\2m? + a) 
(t>l{x) K (floo ^ (51) 

X 

where 0oo and a are some constants. This allows us to assert that there is an 
eigenvalue m\ for which an eigenfunction exists. This is presented in Fig.|Sl to 
some degree of accuracy. 

For the value m* equation H33|l has the following asymptotical behavior: 

(j^'iix) + ^(/)'i « 2A2 {mlf (</)! - ml) (52) 
and the corresponding asymptotical solution is 



(f)i{x) ml + Pi (53) 
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1,0-, 




-1,5-1 



Figure 2: The typical singular solutions f{x) of Eq. The solution is pre- 

sented for i = 4 step. 



where Pi is some constant. 

The next step is finding the first approximation for the fi{x) function. The 
corresponding equation is 

/r + ^/' = /i[0? + Ai(/f-A.?)]- (54) 

The numerical investigation shows that for arbitrary /i there are two differ- 
ent singular solutions which are presented in Fig. [21 Analogously, the singular 
behavior of the solution /i (x) is 

V XiX-Xq 

h{x) « by ^ll>^il. (56) 

\ \l X — Xq 

Evidently, we can suppose that there is the eigenfunction fi{x) for the eigen- 
value fj,i = fil with the following asymptotical behavior: 



n{x) « /oo,i (57) 

where /oo,i is some parameter. 

The next step is the substitution of the first approximations (j)\{x), fl{x) 
into Eqs. (|40|l and H41|l . to find eigenfunctions h\{x),g1{x) corresponding to the 
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eigenvalue E^, 



91 



91 
hi 



El + M - 

-El + M ■ 



ifif + iriY 



(58) 
(59) 



The numerical solutions of these equations are the same as that of the Dirac 
equation for a quantum particle with spin in a potential hole. In our case the 
hole is formed by the gluon fields: Vg// = ct{4''^ + /^)/2- The typical singular 
solutions of these equations are presented in Figs. 13 and 0] 
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E=0.943 
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E=0.93 



Figure 3: The typical singular solutions h{x) of Eqs. (|4U|I and H41|) . 



Further, the eigenfunctions fi{x),hl{x),gl{x) will be substituted into equa- 
tion i^E\i to find the eigenfunction 02 (^) which corresponds to the eigenvalue 
7712 = "^21 ^^'^ SO on. 

The result of these calculations is presented in Figs. EJ |3 |H1 and Table 
ni One can see that there is the convergence: (l)*{x) 4>*{x), fi{x) f*{x), 
h*{x) — > h*{x), g*{x) g*{x), m* m*, /x* fi* and E* ^ E* , where 
f*{x),(p*{x),h*{x),g*{x), are the eigenfunctions and m*,fi*,E* are the eigen- 
values of the nonlinear eigenvalue problem (I25ll - H28|) . 

Finally we would like to mention that the presented solution allows us to 
avoid the Derrick's Theorem jS] which forbids the existence of regular solutions 
for scalar fields. One of the necessary conditions of this Theorem is that the 
solution at infinity tends to a global minimum of the potential energy but in 
our case the solution tends to a local minimum (p* m* , /* ^ (the global 
minimum is (j) — f = fj,). 
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5i 




Figure 4: The typical singular solutions g{x) of Eqs. l|in)l and (|1T|| . 



i 


1 


2 


3 


4 




1.837435109502 


1.61449115 


1.6296365 


1.63265546 




1.4921052897028195 


1.504847 


1.5022397 


1.5030292 


E* 


0.932603482 


0.942702 


0.939357 


0.940497 



Table 1: The iterative parameters m*, fi* and E*. 

5 Properties of the solution 

In this section we would like to describe properties of the derived solution. It is 
easy to see that the asymptotical behavior of the regular solution is 



exp{-a;A/2A2 {m* f} 
4>*{x) « m*+0oo ^- (60) 



exp{-x^(m*)"-Ai ifi*f} 
f [X) « /oo (61) 



eM-J(M-^^Y-iE*f} 
h*{x) « ^— (62) 



eM-xJ(M-^y-{E*f} 
g*{x) « ^ (63) 
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X 



Figure 5: The iterative functions 01,2,3,4(2;). 



-9o 



I M ~a{m*f/2 + E* 



(64) 



M - a {m*y /2 - E* 

where m* , fi*, E* are the eigenvalues derived by the solving Eqs. H25|) - H28|) . 

For determination of the corresponding energy we rewrite the Lagrangian 
l|14(l in the form 

'c) ^ £, + Cnt + (65) 



where Lagrangians for the scalar field interaction Cmt, and spinor field 
are, correspondingly, 



c 



2 A2 

~ T 



4, 



4 

CP 

a (0^0^) (^V) , 



(66) 

(67) 
(68) 



Then the energy density is found as follows: 



sir) =^ 



1 



A. 



y<?!)™0™ (0"<?!)™ - 2<?!)™0g") + (0" </)'') ((^"(Z)") 



(69) 



12 




We will determine the energy density of the quark field from the energy- 
momentum tensor of ip, 



dx°' 



Using ansatz (|2U|) we have 



(70) 
(71) 



Therefore 



G2 



--if* 



iff 



{h* f - {g* f 



{h*f + {g*f 



(72) 

where Ai_2, w* and n* are redefined according to the remark after Eq. (|24ll . and 
the functions (■ • ■ ) are the eigenf unctions. Thus the energy is 



oo 


f 1 










"f - (m*f 



(73) 



-E* 



{h*f + {g*f 



>dx 



47r(/)o/i {\l^2,m*,^l*,E*,a) 
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Figure 7: The iterative functions ft.i,2,3,4(a;)- 



where the quantities r, f*,(t)*,h*,g*, m* , fi*, E*,a are dimensionless according 
to the remark after Eq. (|24|l . The quantity 0(0) ~^ defines the radius of the 
object as a consequence of the asymptotical behavior H60|) - H64|) . The profile of 
the energy density is presented in Fig. O The numerical calculations for the 
integral /i gives 



oo 







(/*)'' + 


\2 

2 








irf- 










a 
~ 2 


(/*)' 




[ih*f - 


(.9* 


f +E* {h*f 





~ 0.395. 
(74) 



Another interesting property of the solution is a total angular momentum 
density Mi. The corresponding operator is 

M, = + s, (75) 

where Li is the operator of the orbital momentum density and is the spin 
operator density, 

si = i"f'^"f^, S2 = i7^7\ S3 = «7^7^ (77) 
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where 7* denote Dirac matrices. The calculation of gives us 

1 

—'I 
2 

{My) 



0. 



(78) 
(79) 



The ^-projection of the full angular momentum is 
M 



1 



-47r / X 



{h*f + {g*r 



dx 



(80) 



47r/2 



where U — 1 and I2 ~ 0.063. Most interesting here is that ijj is not a wave 
function of a quantum particle but it is the quark field which interacts in a non- 
linear way with the scalar fields (which after above-mentioned simplifications 
present gluon fields). Consequently, -0 cannot be normalized to unity, i.e. 



{h^+g')dV^l. 



(81) 



This leads to a very interesting result: the z-projection of the total angular 
momentum of the solution is < This can be related to the spin prob- 

lem of nucleon which raises experimental and theoretical questions regarding 
the contribution of the orbital momentum of the quarks to the total spin of the 
nucleon. One can suppose 0] that the full spin of nucleon has two components: 
one from gluon fields and another from quark field. Our calculations are in 

iNote that for g(0) = 0.2 we have Mz > h/2. 



15 



0,02- 




10 



15 



20 



-0,02- 



X 



^ -0,04- 



-0,06- 



-0,08- 



-0,10- 



Figure 9: The profile of the energy density. 



agreement with this statement: the angular momentum from the quark field is 
non-zero and can be < Ji/2. Unfortunately, for the presented model we cannot 
calculate the term coming from the gluon fields since the gluon fields arc pre- 
sented here by the scalar field which does not have orbital momentum or spin. 
But it can be the aim for future investigations. For the simple model presented 
here, one can try to tune free parameters of the model, (A1.2, M, a, g{0)) in such 
a way so as Mz = h/2, i.e. / (/i^ + 5^) dV = 1, is fulfilled. 

6 Discussion, conclusion and problems 

In this paper, we have presented a model of nucleon based on the non-perturbative 
Heisenberg's quantization technique applied for the SU(3) gauge Yang-Mills 
theory plus quark field. We have shown that under some assumptions and 
simplifications one can reduce the gluon-quark Lagrangian to a scalar- fermion 
Lagrangian. The derived spherically symmetric solution of corresponding equa- 
tions can be considered as a very simplified model of nucleon. Such a "nucleon" 
has a finite mass and spin < h/2. 

The comparison between the "nucleon" and the scalar glueball (which is 
"nucleon" minus quark field) shows us that there is a mass gap which is rclatedyo 
the presence of the energy level E* of the quark field. This value is E* ^ 
because it is the eigenvalue for the Dirac equation. 

Another interesting feature of the "nucleon" is that A4z < h/2. One can 
hope that after some complication of the presented model an additional term 
from the gluon fields will appear and the sum of these two terms gives us the 
spin = fi/2. This could resolve the spin problem of the nucleon. 

Now we would like to list the assumptions and simplifications which are nec- 
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essary for the reduction from the gluon-quark Lagrangian to the scalar-fermion 
Lagrangian: 

• The infinite set of equations connecting Green's functions are reduced 
using the averaged Lagrangian. 

• This reduction is based on the one-function approximation for the 2- and 
4-point Green's functions. 

• It is assumed that the quark field approximately is classical one. 

• It is assumed that the interaction between gluon and quark fields is de- 
termined by their mean-square fluctuations: Cint = (^{^''^ (^^^^^^ 

Also, we would like to mention the following problems arising in the context 
of the considered model: 

• It is necessary to prove the existence of the presented regular solution. 
Mathematically, this means that one should give an exact ground of the 
presented method of the numerical solution of the non-linear eigenvalue 
problem if^ - lf^ . 

• From the physical point of view, it is necessary to generalize this model so 
that it includes the gluon orbital momentum which can give a contribution 
to spins of the glueball and the nucleon. 

• To compare the presented value of the energy with the mass of nucleon 
it is necessary to tune up free parameters of the model. It is not a sim- 
ple procedure since the preliminary numerical analysis shows us that the 
regular solutions exist not for all values of the parameters. 

• It is interesting to investigate the other possible variants of the interaction 
between gluon and quark fields. In this connection, one can mention 
Ref. in which a model of nucleon is derived on the basis of a "chromo 
dielectric model" . One can suppose that there exists such a gluon-quark 
interaction that will lead to the "chromo dielectric model" . 
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